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ABSTRACT 

Although the rapid neutron-capture process, or r-process, is fundamentally important for explaining 
the origin of approximately half of the stable nuclei with A > 60, the astrophysical site of this process 
has not been identified yet. Here we study r-process nucleosynthesis in material that is dynamically 
ejected by tidal and pressure forces during the merging of binary neutron stars (NSs) and within 
milliseconds afterwards. For the first time we make use of relativistic hydrodynamical simulations 
of such events, defining consistently the conditions that determine the nucleosynthesis, i.e., neutron 
enrichment, entropy, early density evolution and thus expansion timescale, and ejecta mass. We find 
that 10 _3 -10 -2 Mq are ejected, which is enough for mergers to be the main source of heavy (A > 140) 
galactic r- nuclei for merger rates of some 10 _5 yr~ 1 . While asymmetric mergers eject 2-3 times more 
mass than symmetric ones, the exact amount depends weakly on whether the NSs have radii of ~15 km 
for a "stiff" nuclear equation of state (EOS) or ~12 km for a "soft" EOS. R-process nucleosynthesis 
during the decompression becomes largely insensitive to the detailed conditions because of efficient 
fission recycling, producing a composition that closely follows the solar r-abundance distribution for 
nuclei with mass numbers A > 140. Estimating the light curve powered by the radioactive decay 
heating of r-process nuclei with an approximative model, we expect high emission in the B-V-R bands 
for 1-2 days with potentially observable longer duration in the case of asymmetric mergers because 
of the larger ejecta mass. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: neu- 
tron 
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1. INTRODUCTION 

The r-process, or rapid neutron-capture process, of 
stellar nucleosynthesis is invoked to explain the pro- 
duction of the stable (and some long-lived radioactive) 
neutron-rich nuclides heavier than iron, which are ob- 
served in stars of different m etallicities as well as i n the 
solar system (for a review, see Arnould et al. |2007 ). De 
spite a growing wealth of observational data (e.g., Sheden 
|et al. |[2008p and although increasingly better r-process 
models with new astrophysical or nuclear physics ingre- 
dients have been developed over decades, the stellar pro- 
duction site(s) of r-process material has (have) not been 
identified yet. All proposed scenarios face serious prob- 
lems. Supernovae (SNe), for example, appear attractive 
because of their potential to explain observational fea- 
tures of the galactic chemical enrichment history (e.g 



Argast et al. 2004). Their nucleosynthesis, however, 
exhibits extreme sensitivity to the detailed conditions in 
the ejecta, whose viability for strong r-processing could 
not b e verified by sophisticated hydrodynamical models 
(e.g ., IHoffman et al. ||2008[ IJanka et al. ||2008[ IRoberts 



et al. 



Wanajo et aL 



Hud epohl et al. ||2010| |Fischer et a 



2011) 



Early in the development of the theory of nucleosynthc- 
sis, an alternative origin of r -process nuclei was proposed 
( Tsuruta fc Cameron |1 965). It relics on the fact that at 
high densities (typically p > 10 10 gem -3 ) matter tends 
to be composed of nuclei lying on the neutron-rich side of 
the valley of nuclear stability as a result of endothermic 



free-electron captures. The astrophysical plausibility of 
this production mechanism as the source of the observed 
r-nuclides has lon g been questioned. It re mained l argely 



unexp lored until Lattimer et al. (1977) and Meyer 



(1989) studied the decompression of cold, neutronized 
matter ejected by tidal effects of a black hole (BH) on 
a neutron star (NS) companion. Recently, special atten- 
tion has been paid to NS mergers because hydrodynamic 
simulations of NS-NS and NS-BH mergers showed that 
a non-negligible amoun t of matter may be ejected (e.g., 
Janka et aL||1999[ [Rosswog et al. ||2004} |Oechslin et al. 
||2007p . 

All pre vious investigations of the ejecta from coalesc- 
ing NSs (JF rciburghaus et al 
Arnould et al. j |2007| |Metzger et al 



Goriely et al. 112005 
SOlOURoBeTte 



et al. 2011 ) were parametrized in one way or another 



which makes their results and conclu sions subject to open 



questi ons. In Goriely et al. ( 2005 ) and Arnould et al. 



(20071 the thermodynamic profiles were constructed by 
a simple decompression model (see Goriely et al. |20TT ), 
but the neutron enrichment (or equivalently the electron 
fraction Y e ) was consistently taken from /3-equilibrium 
assumed to have been achieved at the initial density 
prior to the decompression. It was found that the fi- 
nal composition of the material ejected from the inner 
crust depends on the initial density, at least for the outer 
parts of the inner crust at Pdrip < P < 10 12 g/cm 3 (where 
Pdrip ~4.2xl0 1:L g/cm 3 is the neutron-drip density). For 
the deeper inner-crust layers (p > 10 12 g/cm 3 ), large 



Goriely et al. 



neutron-to-seed ratios drive the nuclear flow into the 
very heavy-mass region, leading to multiple fission re- 
cycling. As a consequence, the resulting abundance dis- 
tribution becomes independent of the initial conditions, 
especially of the initial density. It was found to be in 
close a greement with the solar distribution for A > 140 
nuclei ( |Goriely et aT7~|[2005J |Arnould et al~"|2007 1 . 

Different approaches were t aken to nucleosynthesis cal - 
culations for merger eje cta by Frciburgha us et al. j 1999 1 ; 

Roberts et al. 



effects that become important when the NSs collide. The 
correspon ding NSs have radii of ^15 km (see Bauswein et 
al. |2010). To assess the EOS influence on gross proper- 



Metzger et al. (20101 and 



iMrrin 



their calculations, while the density evolution of the mass 
elements was adopted from hydrodynamical simulations, 
both the initial neutron enrichment and the temperature 
history were considered as free parameters. In particular, 
Y e was chosen in order to obtain, after decompression, an 
r-abundance distribution as close as possib le to the solar 
distribution . This l ed to values of Y f = 0. 1 ( |Freiburghaus| 
|et al. |1999| ) or 0.2 ( |Roberts et al. |2011[ ), corresponding 



to relatively near-surface layers of the inner crust and 
to nuclear flows that are not subject to multiple fission 
cycles. In our simulations, most of the ejecta mass origi- 
nates from the deep layers of the inner crust so that any 
contribution from near-surface layers remains minor. In 
addition, previous studies assumed that the nucleosyn- 
thesis is independent of the initial temperature of the 
ejected material, while we find that the initial temper- 
ature may not only affect the initial composition, but 
potentially also the nucleosynthesis (see Sect. [3]) so that 
special attention should be paid to the detailea temper- 
ature history of the ejected material prior to its free ex- 
pansion. 

The work presented here is based on recent three- 
dimensional relativistic simulations of NS-NS mergers to 
determine the nucleosynthesis-relevant conditions of the 
ejected matter. The hydrodynamical model is described 
in Sect. [2] Section [3] presents the nucleosynthesis re- 
sults for two binaries, a symmetric NS-NS system and an 
asymmetric one. In contrast to previous studies, detailed 
information about the density, Y e , and entropy evolution 
of the ejecta is extracted from the hydrodynamical sim- 
ulations and included in the network calculations. The 
expected electromagnetic emission that is powered by ra- 
dioactive decays following the heavy-element nucleosyn- 
thesis is estimated by a simple, approximative light-curve 
model in Sect.|U Conclusions are drawn in Sect. [5) 

2. HYDRODYNAMICAL MODELS 

Our NS-NS merger simulations were performed with 
a gener al relativistic Smoothed Particle Hydrodyn amics 



scheme ( Oechslin et al. 2007 Bauswein et al. 2010) rep 
resenting the fiuid by a set of particles with constant rest 
mass, whose hydrodynamical properties were evolved 
according to Lagrangian hydrodynamics, keeping Y e of 
fluid elements fixed. The Einstein field equations were 
solved assuming a conformally fiat spatial metric. Bi- 
naries with two mass ratios were modeled, namely sym- 
metric 1.35M0-1.35M0 and asymmetric 1.2Mq-1.5Mq 
systems, both with a resolution of about 550,000 parti- 
cles. The 1.35M -1.35M0 case is of particular interest 
since, according to population synthesis studies and pul- 
sar o bservations, it represent s the most abundant sys- 
tems (Belczynski et al. |2008|). 



ties of the ejecta , we al so performed simulations with the 



For t he results presented here we used the |Shen et ah 
[ (1998) equation of state (EOS), which includes thermal 



Lattimer et al. ( |1991 1 LS220 EOS with an incompress- 
ibility modulus of 220 Me V and a NS radius of ~12km. 
Note that both of our E OSs are in agreement w ith the ob- 
servation of a 2 Mq NS (Demorest et al. 2010 ). Since an 
accurate calculation of the temperature evolution in the 
hydrodynamic simulations is hampered by high initial 
degeneracy and limited numerical resolution, we post- 
processed the temperature of each ejected fluid element 
by assuming an adiabatic flow in the absence of shocks 
and increasing the specific entropy (consistent with the 
shock-jump conditions) when a shock was detected. 

In the case of the symmetric NS binary, about 
3xlO _3 M are found to become gravitationally un- 
bound, whereas about 6xlO _3 M Q are ejected from the 
asymmetric system (using the LS220 EOS, we obtain 
an ejecta mass of ^2xlO~ 3 M0 and 6xlO _3 M Q respec- 
tively). The ejected "particles" (i.e., mass elements) orig- 
inate mostly from two different regions in the inner crust 
of the initial stars. For the symmetric model ^75% of 
the material are squeezed out from the contact interface 
of the NSs. The remaining 25% are ejected from the 
near-surface regions close to the orbital plane, see inset 
of Fig. [II where the arrows indicate by trend which ini- 
tial densities the majority of fluid elements corresponds 
to. Note that the "contact" interface of the asymmetric 
merger is different from the symmetric case. While in the 
latter the two NSs collide in a shearing way and violently 
along a broad area, the impact of the tidally stretched 
lower-mass companion in the asymmetric case happens 
in a more grazing way such that a conical "nose" of 
the lower-mass star digs off matter from the higher-mass 
object. For this reason, in the asymmetric system the 
contact region preferentially contributes to ejecta with 
low initial densities. As seen in Fig. [T] all ejecta have 
low Y e ranging between 0.015 and 0.050 (about 30% of 
the mass have Y e ~ 0.015). In the asymmetric case a 
small amount of ejected material starts out from densi- 
ties above 2xl0 14 g/cm 3 and has Y e - values up to 0.07. 
For the first 27 ms the density history is consistently fol- 
lowed by the numerical simulation. Afterwards the es- 
caping ejecta are assumed to expand freely with constant 
velocity. The radii of the ejecta clumps thus grow lin- 
early with time t and consequently their densities drop 
like 1/t 3 . Note that due to a lack of resolution the dy- 
namics and mass of unbound crust material with initial 
densities p < 10 13 g/cm 3 cannot be reliably calculated 
in the hydrodynamic model. However, compared to the 
1O _3 -1O _2 M0 of inner crust matter its contribution to 
the ejecta remains small. 

The ejected matter is initially cold, but most of it gets 
shock-heated during the ejection to temperatures above 
1 MeV. Its composition is then determined by nuclear 
statistical equilibrium (NSE). When the drip density is 
reached during expansion, most of the matter has cooled 
below 1 MeV and the NSE composition has frozen out. 
As soon as the temperature has dropped below 10 10 K, 
further changes of the composition are followed by a full 
network calculation (as detailed below) , and the temper- 
ature evolution is determined on the basis of the laws 
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Fig. 1. — Histograms of fractional mass distribution of the ejecta for the 1.35-1.35Mq NS merger (upper row) and the 1.2-1.5M0 binary 
(lower row) as functions of density p (relative to the saturation density ps — 2.6x 10 14 g/cm 3 ; left) and of electron fraction Y e (middle) that 
the ejected matter had at its initial NS location prior to merging. The right panels show the fractional mass distributions as functions of 
the final entropy S per nucleon when the matter starts its free expansion. In the insert on the left panels the dots mark positions of mass 
elements that get ejected later. The locations are given in projection on the orbital plane at the time when the stellar collision begins. 



of thermodynamics, allowing for possible nuclear heat- 
in g through ^-de cays, fission, and a-decays, as described 



Meyer (1989 1 . 



3. NUCLEOSYNTHESIS 

The reaction network includes all 5000 species from 
protons up to Z=110 lying between the valley of j3- 
stability and the neutron-drip line. All fusion reactions 
on light elements that play a role when the NSE freezes 
out are included in addition to radiative neutron cap- 
tures and photodisintegrations. The reaction rates on 
light species are taken from the NETGEN library, which 
includes all the latest co mpilations of expe rimentally de- 



termined reaction rates (Xu et al. 2011). Experimen- 



tally unknown reactions are estimated with the TALYS 



2008[) on the basis o f the HFB- 
1 Goriely et ^[2010] ). On top of 



code (Goriely et al. 
21 nuclear mass mode 

these reactions, fission and /^-decays are also included, i.e 
neutron-induced fission, spontaneous fission, ^-delayed 
fission, photofission, as well as /3-delayed neutron emis- 
sion. The /3-decay processes ar e taken from the updated 



version of the Gross Theory (Tachibana et al. 19901 



based on the HFB-14 Q- values, whereas all fission pro 
cesses are estimated on the basis of the HFB-14 fission 
path and the full c alculation of the corresponding barrier 
penetration (Goriely et al. 20091. The main fission re- 



gion is illustrated inFig.[2[ The fission fragment distribu- 
tion is taken from Kodoma & Takahashi ( 1975 1, and the 
fragment mass- and charge-asymmetry are derived from 
the HFB-14 prediction of the left- right asymmetry at the 
outer saddle point. Due to the specific initial conditions 
of high neutron densities (typically N n ~ 10 33_35 cm -3 
at the drip density), the nuclear flow during most of the 
neutron irradiation will follow the neutron-drip line. For 
these nuclei at T > 2-3 xlO 9 K, (n,2n) and (2n,n) reac- 
tions are faster than (7,11) and (11,7) reactions and must 
be included in the reaction network. The (n,2n) rates 
are estimated with the TALYS code and the reverse rates 
from detailed balance expressions. 

For drip-line nuclei with Z > 103, fission becomes ef- 
ficient (Fig. I2J) and recycling takes place two to three 
times before the neutrons are exhausted, as shown in 
Fig. ^ by the time evolution of the mass number (A) 
mass-averaged over all the ejecta. After several hundred 
ms, when neutrons get exhausted by captures (N n ~ 
10 20 cm -3 ), n-captures and /3-decays compete on similar 
timescales and fashion the final abundance pattern be- 
fore the nuclear flow becomes dominated by /3-decays (as 
well as fission and a-decays for the heaviest species) back 
to the stability line. The average temperature remains 
rather low during the late neutron irradiation, around 
0.5 GK (Fig. [3]) so that photoreactions do not play a ma- 
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Fig. 2. — Representation of dominant fission regions in the (TV, Z) 
plane. Nuclei for which spontaneous fission is estimated to be faster 
than /3-decays are shown by full squares, those for which /3-delayed 
fission is faster than /3-decays by open squares, and those for which 
neutron-induced fission is faster than radiative neutron capture at 
T = 10 9 K by diamonds. 

jor role. 

The final mass-integrated ejecta composition is shown 
in Fig. [4J The A = 195 abundance peak related to 
the N = 126 shell closure is produced in solar distri- 
bution and found to be almost insensitive to all input 
parameters such as the initial abundances, the expansion 
timescales, and the adopted nuclear models. In contrast, 
the peak around A = 140 originates exclusively from 
the fission recycling, which takes place in the A ~ 280- 
290 region at the time all neutrons have been captured. 
These nuclei are predicted to fission symmetrically as vis- 
ible in Fig. [4] by the A ~ 140 peak corresponding to the 
mass-symmetric fragment distribution. It is emphasized 
that significant uncertainties still affect the prediction of 
fission probabilities and fragment distributions so that 
the exact strength and location of the A ~ 140 fission 
peak (as well as the possible A = 165 bump observed 
in the solar distribution) depend on the adopted nuclear 
model. 

While most of the matter trajectories are subject to a 
density and temperature history leading to the nuclear 
flow and abundance distribution described above, some 
mass elements can be shock-heated at relatively low den- 
sities. Typically at p > 10 10 g/cm 3 the Coulomb effects 
shift the NS E abundance distrib ution towards the high- 
mass region ( Goriely et al. |2011 1, but at lower densities, 
the high temperatures lead to the photodissociation of all 
the medium-mass seed nuclei into neutrons and protons. 
Nucleon recombination may occur during the decompres- 
sion provided the expansion timescale of the trajectories 
is long enough. For a non- negligible amount of ejected 
material, this recombination is indeed inefficient so that 
light species (including D and 4 He) are also found in the 
ejecta (Fig. |4| . The final yields of A < 140 nuclei re- 
main, however, small and are not expected to contribute 
to any significant enrichment of the interstellar medium 
compared to the heavier r-elements. 

4. ELECTROMAGNETIC COUNTERPARTS 

Radioactive power through /3-decays, fission processes 
as well as late-time a-decays will heat the expanding 
ejecta and make them radi ate as a "macro-nova" (JKulka- 



2005) or "kilo-nova" (Metzger et al. 2010I associ 



ated with th e ejection of nuc leosynthesis products from 



the merger ( Li fc Paczynski~| [T998 ) . The time evolu- 
tion of the corresponding total mass-averaged energy re- 
lease rate available for heating the ejecta (i.e., energy 
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Fig. 3. — Time evolution of the total radioactive heating rate 
per unit mass, (Q), mass number (A), and temperature (T) (all 
mass-averaged over the ejecta) for the 1.35-1.35M0 (solid lines) 
and 1.2-1.5M0 (dotted lines) NS mergers. 
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Fig. 4. — Final nuclear abundance distributions of the ejecta from 
1.35— 1.35M0 (squares) and 1.2-1.5Mq (diamonds) NS mergers as 
functions of atomic mass. The distributions are normalized to the 
solar r-abundance distribution (dotted circles). 

escaping in neutrinos is not considered) is plotted in 
Fig. [3] for both the 1.35-1.35M and 1.2-1.5M Q bina- 
ries. \Vhile (Q(t)) and the average temperature evolu- 
tion differ only slightly between both NS-NS systems, the 
ejecta masses M e j and mass-averaged expansion veloci- 
ties v exp differ considerably. While we find for the sym- 



metric system v exp ~ 0.31c (c being the speed of light) 
and M C j m 3 x 1O~ 3 M , corresponding to a total heat- 
ing energy of -E^cat ~ 2 x 10 49 erg or 3.4MeV/nucleon, 
the numbers for the asymmetric case are v cxp s=s 0.23c, 
M e j w 6 x 10" 3 M Q and E hcat w 4 x 10 49 erg (again 
3.4MeV/nucleonjM This must be expected to lead to 
significant differences in the brightness evolution of the 
kilo-nova because its peak bolometric luminosity scales 

with I/p ea k oc VcxpM- and, for free expansion (u exp = 
const), is reached on a time-scale i poa k oc v CX p M / 



( |Metzger et al. ||2010[ |Arnett |1982[ ). 

We calculate an approximation of the light-curves of 



In the simulations with the LS220 EOS we obtain v c 



0.28c 



for the symmetric and v B 



i 0.24c for the asymmetric binary. 
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Fig. 5. — Photon luminosities of the expanding NS merger ejecta 
caused by radioactive decay heating for the 1.35-1.35M0 (solid 
lines) and 1.2— 1.5Mq (dashed lines) binaries. The upper, long- 
duration lines are the bolometric luminosities, the sequences of 
short-duration peaks correspond to the emission in the blue, visual, 
and red wavebands (at wavelengths of 445, 551, 658 nm; from left 
to right). 

such events by employing; a semi -analytic, simplified one- 



zone model of Kulkarni 
Li fc Paczyhski ||liW8) ; 



(2005) (see also Arnett 1982 



assumi ng that asymmet ries o 



the emission remain modest (cf. Roberts et al. 2011) 



The EOS includes photons as well as nuclei and elec- 
trons, whose recombination fraction with decr easing tem- 
perature is appro ximated by that of 56 Ni (Arnaud & 
Rothenflug ||1985[ ). The opacity k = 0AZ/Acm z g~ l for 
Thomson scattering is used, taking mass-averaged Z and 
A from the nucleosynthesis yields (Fig. pi) and ignoring 
electro n recombination as suggested by |Roberts et al. \ 

poiTb. 



Our results for the bolometric light curves L^ \(i) and 
the B-, V-, and R-band luminosities vL v {t) are displayed 
in Fig.[5j While we expect significant emission in the cho- 
sen wavelength bands up to (^£j,) pe ak ~ 4 x 10 41 ergs -1 
for about one day in the case of a 1.35-1.35Mq merger, 
the 1.2— 1.5Mq system produces sizable BVU-radiation 
nearly twice as long. 

5. CONCLUSIONS 

Using relativistic NS merger models to determine the 
nucleosynthesis-relevant conditions self-consistently, we 
confirm that decompressed NS matter ejected dynami- 
cally during the stellar collision and shortly afterwards 
is an extremely promising site for robust, strong r- 
processing. Matter from the inner crust of the coalescing 
NSs, which dominates the ejecta by far, produces an r- 
abundance distribution very similar to the solar one for 



tribution c ould originate from the outer crust ( Goriely et 



al. 2011[), but too little of such matter gets ejected to 



explain the solar proportion of light-to-heavy r-process 
material. However, significant amounts of A < 140 nuclei 
might be produced in the outflow of a BH-torus system 
formed after the NS merger (Wanajo & Janka 2011). 

The underlying nuclear mechanisms differ significantly 
from those at action in SN scenarios. In particular fission 
plays a major role in recycling heavy material. The simi- 
larity between predicted and solar abundance patterns as 
well as the robustness of the prediction against variations 
of input pa rameters, which we have shown in|Goriely et| 



al. ( 2005 1 , have demonstrated here in comparison of 



symmetric and asymmetric NS-NS mergers, and will fur- 
ther elaborate on in a forthcoming paper, make this site 
one of the most promising, deserving further exploration 
concerning hydrodynamics, galactic chemical evolution, 
nucleosynthesis, nuclear physics, and astronomical con- 
sequences. 

Fully relativistic simulations including neutrino trans- 
port and magnetic fields, with good resolution of the 
inner (and outer) crust layers of the merging NSs as- 
sured by adaptive refinement, are needed to corrobo- 
rate the ejecta conditions found in our work. With 
our yield of ^(3-5) xlO -5 M© per event of 151 - 153 Eu, a 
nearly-pure r-process element, the origin of all galactic 
Eu from NS-NS mergers would require an event rate of 

2-3)xl0 _5 yr _1 , fully compati ble with pre sent best es- 



timates (see, e.g., |Belczynski et al" 



also mass loss during the evolution o 



2008]). However, 
the merger rem- 



nant from a hypermassive NS to a massively accreting 
BH-torus system as well as NS-BH mergers with higher 
ejecta masses might add to the production (e.g., Wanajo 
& Janka 2011, Caballero et al. 2011) and deserve more 
detailed investigations. From the nuclear point of view, 
this site also implies new challenges because it involves 
the formation of neutron-drip nuclei, for which /3-decay, 
neutron capture, and fission rates need to be determined. 
Astronomically, the discov ery of electromagnetic radia- 
tion from kilo- nova events (iRoberts et al. 120111 iMetzg 



2010) could mean the first in-situ observation 



et al. 

freshly produced r-process material. 
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